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Summary  and  Conclusions 


Motivation 

•  Can  Nanotubes  conduct  heat  well? 


•  Which  Phonon  Modes  limit  Thermal  Conduction? 


•  Can  interacting  Nanotubes  Spin  and  Twist? 
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•  Is  nanostructured  carbon  foam  the  ultimate  structural 
and  thermal  material? 


Morphology  of  carbon  nanotubes 


Structure  of  carbon  nanotubes: 
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multi-wall  tube 
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rope/bundle  of  tubes 


•  Classification  of  carbon  nanotubes: 

zigzag 


*:metal  •  tsemiconductor 


armchair' 


R  =  mo!  +  nao 
defines  uniquely  the 
(m,n)  nanotube 


•  Types  of  achiral  nanotubes: 
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Computational  Techniques 


Total  Energy  Calculations 

-  A  b  initio  Density  Functional  Formalism 

*  Plane- wave  code  with  energy  cutoff  of  50  Rydbergs 

*  Soft  nonlocal  pseudopotentials  and  Ceperley- Alder  exchange- 
correlation  potential 

*  200,000  plane  waves  in  the  basis  set 

—  Parametrized  Linear  Combination  of  Atomic  Orbitals  (LCAO) 
Formalism 

*  Parametrization  based  on  ab  initio  Density  Functional  re¬ 
sults 

*  Very  fine  k-points  mesh  for  the  inter-tube  interaction 

*  O(N)  technique  for  very  large  systems 

*  Ideally  suitable  for  massively  parallel  computation 

—  Tersoff  Potential 

Thermodynamics:  Monte  Carlo  simulations 

Time-dependent  processes: 

Molecular  dynamics  simulations 

Transport:  Landauer-Biittiker  Formalism 

-  Conductance  G(E)  evaluated  by  using 

G(E)  =  ^T(E)  =  ^T{E) 


Thermal  Conductivity  in  Nanotubes 
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•  Experimental  Observation  for  a  Nanotube  “Mat” 


J.  Hone,  M.  Whitney,  A.  Zettl, 

Synthetic  Metals  103,  2498  (1999) 

—  Raw  data  (T  —  300  K):  Ara0.7  W/m-K 

-  Extrapolated  data  (T  =  300  K):  Aral,  800  —  6, 000  W/m-K 


h-ajt-w* 


Approaches  and  Pitfalls 


Direct  Molecular  Dynamics  Approach 
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—  Temperature  profile  in  a  (10, 10)  nanotube 


z  [A] 


•  Pitfalls  of  Direct  Molecular  Dynamics: 

—  L  «  /«100  nm 

=>■  velocity  rescaling  mimics  “grain  boundary” ,  limits  l 

—  Large  statistical  uncertainty 

—  T(z)  not  linear 

—  dT /dz  »  0  =>■  spatial  variation  of  thermal  conductivity 


Equilibrium  Molecular  Dynamics  Approach 

-  Green-Kubo  formula  for  the  Navier-Stokes  thermal  conduc¬ 
tivity: 

x=mjr1?<3w-w>* 

J(*)  =  4  E  =  E  Vi  A  ei  -EE  Tijfaj  ■ 

®  i  i  i  j(^i) 

—  Heat  flux  autocorreltation  function  in  a  (10, 10)  nanotube 


Pitfalls  of  Equilibrium  Molecular  Dynamics: 

—  Slow  convergence 

—  Extensive  ensemble  averaging  crucial 

—  Autocorrelation  function  requires  large  data  storage 


Non-Equilibrium  Molecular  Dynamics  Approach 

A  =  lim  A(Fe)  ,  A(Fe)=  lim  — ~^FeT’5  > 

Fe->0  V  e;  ’  y  eJ  t-> oo  FeTV 

Fe  —  thermal  force,  a  —  Nose-Hoover  thermostat  multiplier 


AFj  =  AejFe  —  E  %(r^--Fe)  +  —  E  E  tjk(rjk-Fe)  - 
j(&)  ^  i  H&) 


Axial  heat  flux  in  a  (10, 10)  nanotube  for  Fe  —  0.2  A  1 
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Heat  transport  dependence  on  Fe 


Fe  [A-1] 


Pitfalls  of  Non- Equilibrium  Molecular  Dynamics: 

—  Fe  — >  0  limit  requires  careful  evaluation 

—  Careful  coupling  to  thermostat  required 


Temperature  Dependent  Thermal  Conductivity 


•  Thermal  conductivity  of  an  isolated 
(10,10)  carbon  nanotube 


S 

| 


4xl04 
3xl04 
2xl04 
lxlO4 
OxlO4 

0  100  200  300  400 
T[  K] 


•  Thermal  conductivity  of  various  carbon  allotropes 
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Effect  of  Inter- Wall  Interactions  on  Thermal  Conductivity 


► 

•  Thermal  conductivity  of  a 
(5,5)@(10,10)  double-wall  carbon  nanotube 


•  Thermal  conductivity  of  bundled  (10,10)  nanotubes 


►  Effect  of  Isotope  Mixture  on  Thermal  Conductivity 
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Low-Frequency  Phonon  Modes 
of  Nanotubes 


•  Equilibrium  structure  of  bundled  rigid  (10,10)  tubes 
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►  Equilibrium  Structure  of  Bundled  and  Multi- Wall  Tubes 
•  Interaction  between  neighboring  (10,10)  tubes 


—  Two-dimensional  oblique  lattice. 


•  (5,5)@(10,10)  Double- Wall  Tube  Geometry 


[. Electronic  and  structural  properties  of  multiwall  carbon  nan¬ 
otubes ,  Young-Kyun  Kwon  and  David  Tomanek,  Phys.  Rev.  B 
58,  R16001  (1998)] 


(5,5)@(10,10)  Double- Wall  Tube 

Ar  =  0  A  z  =  0 

\ Pin  0  Pout  —  9 


—  Dependence  of  Energy  on  Off-Axis  Displacement  Ar 
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►  Libration  and  Rotation  Dynamics 
•  (10,10)  tube  “ropes”  [PRB  58,  R1SS14  (1998)} 


—  Librational  motion  at  T  <  Tom 

—  Libration  frequency  z/«12  —  60  cm-1 

—  Experimental  observation  (W.  Holmes  et  al.)\ 

*  Infrared  modes  at  22, 40  cm-1 

*  These  modes  disappear  at  Tm 30  —  180  K 


•  (5,5)@(10,10)  double- wall  tube  [PRB  58,  R16001  (1998)} 


—  Librational  motion  at  T  <  Tom 

—  Libration  frequency  u  ~  33  cm" 
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►  Soft  Phonon  Modes  of  Single- Wall  Tubes 


«■ 


L  =  L  o-AL 


straight 
pinch  mode 
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►  Thermal  Expansion  of  Single- Wall  Tubes 


P  Thermal  contraction  dominates  over  bond  expansion 


at  low  temperatures 


►  Thermodynamics  of  Orientational  Melting 
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•  Floppy  tubes  -  unrealistic  model 
(Torsionally  decoupled  “layers”  within  tube) 


—  Orientational  melting  at  Tom~55  K 
—  Librational  motion  for  TfjToM 
—  Free  rotation  at  T Tom 

•  Rigid  tubes  -  realistic  model 
(torsional  coupling  between  “layers”  within  tube) 

—  Torsion  energy 

AEtorM(Aip/AL)^K  (ff)2 
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*  Torsion  is  thermally  activated  duriner  synthesis  -fwls+tn  fc-WP-iM* 


Orientational  melting 
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T=  cryr0H) 

*  Orientational  disclo cations,  introduced  during  synthesis,  are 
pinned  for  T^Tom- 

*  Axial  motion  of  orientational  dislocations  corresponds  to 
tube  rotation. 

*  Orientational  melting,  described  by  free  “diffusion”  of  ori¬ 
entational  dislocations  along  tubes,  occurs  for  T^Tom~100  K. 

*  Orientational  dislocations  (twistons)  are  scattering  centers 
for  electrons  — »  strong  effect  on  transport  (ID  system!). 
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Morphology  of  Hybrid  Graphitic  Foam 


•Relationship  to  graphite 


refold... 


New  Hybrid  in-between 
Diamond  and  Graphite 


•Folding  Process 


64  carbon  atoms 


unit  cell 


•Hybrid  Graphitic  Foam 


Terrace  size  —>  oo:  graphite 
Terrace  size  0:  diamond 

Present  study: 

• Optimized  structure 
•Elastic  properties 
•Electronic  structure 


sp2  carbon  atoms 


carbon  atoms 


Optimized  Structure 


Elastic  Properties 


New  graphitic  foam  Graphite  Diamond 

P  2.48  g/cm3  2.27  g/cm3  3.51  g/cm3 

C-j  i  (along  a  axis)  9.02  Mbar  12.3Mbar  11.29Mbar 

C22  (along  b  axis)  9.42  Mbar 

C33  (along  c  axis)  0.81  6  Mbar  0.34  Mbar 

Bulk  modulus  0.792  Mbar  0.326  Mbar  4.69  Mbar 


Electronic  Structure 


Brillouin 

Zone 
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Graphitic  foam  is  metallic:  N(EF)=1.0  states/eV/spin 


Density 
of  States 


0  5  10  15 

states/eV  spin 


Summary  and  Conclusions 


We  studied  the  structural  rigidity,  phonon  modes,  and  thermal 
conductivity  of  carbon  nanostructures:  nanotubes  and  graphitic 
foam. 

We  found  that 

-  The  unusually  high  thermal  conductance  of  nanotubes 
results  from  the  large  phonon  mean  free  path. 

-  Weak  inter-tube  coupling  results  in  a  very  soft  librational 
motion. 

-  Orientational  dislocations  are  frozen  in  during  the  syn¬ 
thesis  of  interacting  nanotubes.  Relatively  free  diffusion  mo¬ 
tion  of  orientational  dislocations  occurs  for  T^Tqm  ~160  K, 
marking  the  onset  of  orientational  melting. 

-  Hybrid  graphitic  foam  offers  an  unsusual  combina¬ 
tion  of  properties:  low  mass  density,  high  structural  rigid¬ 
ity,  possibly  high  electric  and  thermal  conductivity. 


